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Hybrid baryon states are described in quark potential models as having explicit excitation of the 
gluon degrees of freedom. Such states are described in a model motivated by the strong coupling 
limit of Hamiltonian lattice gauge theory, where three flux tubes meeting at a junction play the 
role of the glue. The adiabatic approximation for the quark motion is used, and the flux tubes 
and junction are modeled by beads which are attracted to each other and the quarks by a linear 
potential, and vibrate in various string modes. Quantum numbers and estimates of the energies of 
the lightest hybrid baryons are provided. 
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Experiments at new electron-scattering facilities, such 
as those in Hall B at TJNAF, are expected to produce 
many new baryon states, including those which are de- 
scribed in quark potential models as having explicit exci- 
tation of the gluon degrees of freedom. Low- lying baryon 
states present in analyses of ttN elastic and inelastic scat- 
tering, such as the Pn Roper resonance A^(1440) which 
has the same quantum numbers as the nucleon, have been 
proposed j^J^ as hybrid candidates. This is based on 
extensions of the MIT bag model 1^ to states where a 
constituent gluon in the lowest energy transverse electric 
mode combines with three quarks in a color octet state 
to form a colorless state, and on a calculation using QCD 
sum rules . Hybrid baryons have also been constructed 
recently in the large- A'^c limit of QCD Q. 

With the assumption that the quarks are in an S-wave 
spatial ground state, and considering the mixed exchange 
symmetry of octet color wavefunctions of the quarks, 
bag- model constructions show that adding a — 1+ 
gluon to three light quarks with total quark-spin 1/2 
yields both N {I = ^) and A (/ = |) hybrids with 

— i^, 1^. Quark-spin 3/2 hybrids are A^ states with 

= i^, 1"'^, and 1^. Energies are estimated using 
the usual bag Hamiltonian plus gluon kinetic energy, ad- 
ditional color-Coulomb energy, and one-gluon exchange 
plus gluon-Compton 0(as) corrections. Mixings between 
g"^ and q'^g states from gluon radiation are evaluated. 
If the gluon self-energy is included, the lightest A^ hy- 
brid state has = \^ ^^i^ ^ mass between that of the 
Roper resonance and the next observed = state, 

A^(1710). A second = i+ A^ hybrid and a J-^" = |^ 
A^ hybrid are expected to be 250 MeV heavier, with the 
A hybrid states heavier still. A similar mass estimate of 
about 1500 MeV for the lightest hybrid is attained in the 
QCD sum rules calculation of Ref. 

For this reason there has been considerable interest in 
the presence or absence of light hybrid states in the Pn 
and other positive-parity partial waves in -kN scatter- 



ing. Interestingly, quark potential models which assume 
a q'^ structure for the Roper resonance predict an en- 
ergy which is roughly 100 MeV too high, and the same is 
true of the A(1600), the lightest radial recurrence of the 
ground state = | A(1232). Furthermore, models of 
the electromagnetic couplings of baryons have difficulty 
accommodating the substantial Roper resonance photo- 
coupling extracted from pion photoproduction data 
Evidence for two resonances near 1440 MeV in the Pn 
partial wave in tt A^ scattering was cited , which would 
indicate the presence of more states in this energy re- 
gion than required by the q'^ model, but this has been 
interpreted as due to complications in the structure of 
the Pii partial wave in this region, and not an additional 
physical state ||] . 

This letter will show that when the glue in a baryon is 
given the flux-tube structure expected |9j from an expan- 
sion around the strong-coupling limit of the Hamiltonian 
formulation of lattice QCD (HLGT), that the lightest 
hybrid baryons have similar good quantum numbers, but 
substantially higher energies and different internal struc- 
ture than predicted using bag models. This structure 
of the glue, where the gluon degrees of freedom collec- 
tively condense into flux-tubes, is very different from the 
constituent-gluon picture of the bag model and large- A^c 
constructions. The basis of this model is the assump- 
tion that the dynamics relevant to the structure of hy- 
brids is that of confinement. At large interquark separa- 
tions, there is evidence that the flux-tube model hybrid 
potential is consistent with that evaluated from lattice 
QCD l^-jlltl, whereas the adiabatic bag model does not 
reproduce the lattice results there ||l^. It has also been 
shown |l^ that a constituent-gluon model is not able to 
reproduce these lattice results . 

The model used here to describe the glue is the non- 



relativistic flux-tube model of Ref. [^14|, coupled with 
the adiabatic approximation, where the quarks do not 
move in response to the motion of the glue (apart from 
moving as a rigid body in order to maintain the center-of- 
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mass position). Although exact only in the heavy-quark 
limit, this approximation has been shown to be good for 
light-quark mesons in the flux-tube model ||l^. In the 
strong coupling limit of HLGT, flux lines (strings) with 
energy proportional to their length play the role of the 
glue, which are modeled here by equal mass beads with 
a linear potential between nearest neighbors. The total 
mass of all of the beads is given by the energy in the flux 
lines, which is fixed by the string tension. 

Perturbations to the strong-coupling limit are provided 
by the plaquette operator, which moves the flux lines, 
and so the beads, perpendicular to their rest positions. 
Global color gauge invariance requires that the three flux 
lines which emanate from the quarks must meet at a 
junction, which is also modeled by a bead. Lattice QCD 
studies pq] support confinement from a Y-shaped linear 
string rather than pairwise linear strings. As a single 
plaquette operator cannot move the junction and leave 
the Y-string in its ground state, the junction bead is given 
a mass higher than those on the strings. 

The ground state energy of this configuration of beads 
representing the Y-string for definite quark positions 
Ti defines an adiabatic potential VB(ri,r2,r3) for the 
quarks, which consists of the string energy bj^i h, where 
b is the string tension and U is the magnitude of the vector 
\i from the equilibrium junction position to the position 
of quark i, plus the zero-point energy of the beads. The 
energy of the first excited state defines a new adiabatic 
potential Vff (ri, r2, ra). The simplest model incorporat- 
ing the essential degrees of freedom has one bead per 
string, plus a junction bead. The mass of a bead on 
string i is taken to be proportional to bU, and the junc- 
tion bead is given a higher mass. This model has nine 
degrees of freedom: three string bead transverse posi- 
tions (^i) within, and three (zi) out of the plane of the 
quarks, and three junction positions, in (x, y) and out 
(z) of this plane. 

A nonrelativistic Hamiltonian is constructed in the 
small oscillations approximation, and has the form 
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Tj(r) + Tb(e.,i.0 + 71.b- 
+V;{v) + V,,{U,z.,), 
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(0.1) 



where in the potential the string-bead (b) and junction- 
bead (j) coordinates decouple, but there are terms Tbb 
and Tbj in the kinetic energy which couple these mo- 
tions. This Hamiltonian is corrected for the center-of- 
mass (cm.) motion due to the string motion, by allowing 
the quarks to move rigidly to maintain the cm. position. 
This gives effective masses similar to reduced masses to 
the string and junction beads, which depend on the quark 
masses. 

Diagonalization of the resulting coupled-oscillator 
Hamiltonian for a wide variety of quark configurations 
shows that neglect of the coupling terms Tbb and Tbj 
does not significantly affect the lowest string energies. 
Examination of the eigenfunction for the lowest energy 
mode shows that it is always predominantly in-plane 



TABLE L Lowest three non-interacting (NI) and ex- 
act frequencies in GeV for three quark configurations, for 
rrij = rriquark = 0.33 GeV, and 6 = 0.18 GeV^. 



;,(fm) 


£i(NI) 


El 


S2(NI) 


E2 


-E3(NI) 


Ea 


0.5, 0.5, 0.5 


0.614 


0.607 


0.614 


0.607 


0.869 


0.828 


0.5, 0.5, 0.1 


0.623 


0.616 


1.069 


0.985 


1.069 


1.005 


0.5, 1.0, 0.1 


0.520 


0.483 


0.544 


0.534 


0.544 


0.590 



junction motion. This is illustrated for three quark posi- 
tions (specified by the lengths li of the three strings) in 
Table! 

Accordingly, the string ground state and first excited 
state adiabatic surfaces may be found by allowing the 
junction to move, while the strings connecting the junc- 
tion to the quarks follow without excitation. This gives 
the junction an effective mass which in the limit of a large 
number of beads becomes 



M^s = bJ2h 



4^,(W, + M,) 



(0.2) 



where the second term is the center of mass correction 
with quark masses Mi . The potential is the string tension 
times the length of the lines connecting the displaced 
junction to the quarks, so that the junction Hamiltonian 
is 
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influx = ^M^ar^ 



(0.3) 



The adiabatic surfaces are found numerically via a varia- 
tional calculation. This is made necessary by a singular- 
ity in the small oscillations expansion for quark configu- 
rations where the triangle made by joining them contains 
an angle 9i of at least 120°, so that the equilibrium po- 
sition of the junction is such that the corresponding dis- 
tance k between the junction and the quark is zero. This 
variational calculation agrees with the analytic small os- 
cillations results when the li are all large, but shows 
the small oscillations approximation to be poor for some 
other configurations. It was shown in Ref. 10 that when 
the small oscillations approximation is removed, as has 
been done here, hybrid meson masses go down, and when 
the adiabatic approximation is removed, as has been done 
partially here, hybrid meson masses go up. This same 
behavior is seen here. 



As Tj(r) -I- Vj{r) is even under 



-z, one of the 



three first excited modes of the junction always involves 
motion along fj^ = z. Analytic results in the small oscil- 
lations approximation show that the frequency for mo- 
tion along fj^ is always higher than those in the plane of 
the quarks. Trial wavefunctions for the ground and first 
excited states are taken to be the anisotropic harmonic 
oscillator wavefunctions 
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exp {- [{a+i)^ ■ r)2 + {a^fj_ ■ + (a^z)^] /2} 
^H{r)^V2a^fj_-r-^B{r), (0.4) 

with variational parameters a_, a+, a^, and 0, the lat- 
ter giving the direction fj_ in the plane of the quarks 
of the lowest-energy oscillation relative to a body-fixed 
axis in that plane. For every configuration of the quarks, 
specified by the magnitudes of the Jacobi coordinates p, 
A, and the cosine of the angle 9px between them, the 
ground (Vg) and first excited state (Vh) string energies 
are independently minimized. 

In order to compare to the relativized model calcula- 
tion of baryon masses of Ref . , hybrid baryon masses 
are found by allowing the quarks to move in a confining 
potential given by the linear potential b J2i h of Ref. , 
plus Vh — Vb- The Coulomb potential from one-gluon ex- 
change is assumed the same in conventional and hybrid 
baryons, and spin-dependent terms are neglected. Quark 
wavefunctions are expanded in a large oscillator basis, 
and energies for baryons (linear confinement) and hy- 
brids (linear plus Vh~Vb confinement) composed of light 
quarks are evaluated. When added to the spin-averaged 
mass of the TV and A which is 1085 MeV, hybrids with 
quark orbital angular momenta Lq = 0, 1, 2 have masses 
1980, 2340 and 2620 MeV respectively. Hyperfine (con- 
tact plus tensor) interactions split the N hybrids down 
and the A hybrids up by similar amounts, so that the 
N hybrid mass becomes 1870 MeV. The model error on 
this mass is estimated to be less than ±100 MeV. This 
lightest {Lq = 0) hybrid level is substantially higher than 
the roughly 1.5 GeV estimated from bag model and QCD 
sum rules calculations. 

This calculation also shows that the size of the quark 
core of hybrid baryons is roughly 20% larger for hybrid 
baryons than for conventional baryons. 

The hyperfine interaction derived from the Coulomb 
interaction between the quarks has the same sign in con- 
ventional and hybrid baryons, so that the hybrids with 
quark-spin- 1 will be heavier than those with quark-spin- 
i. Given the increased size of the hybrid states, the 
quark-spin- i hybrids are not expected to split as far from 
the spin averaged level as in the usual baryons. 

For every set of quark positions the potential in 
which the junction moves is anisotropic, which means 
that the solutions of the Schrodinger equation for the 
junction motion do not have definite angular momentum. 
However, in the absence of the adiabatic approximation 
the combined wavefunction of the quark and junction 
motions must be a state of good angular momentum. 
Although it is possible to project states of good angu- 
lar momentum out of the combined junction and quark 
motion states, this is technically difficult and will not 
be reported on here. Instead, an intuitive argument is 
given to justify the expectation that the total orbital an- 
gular momentum of the lowest-lying hybrid baryons (H) 



is unity. 

The H hybrid wavefunction in Eq. ( 3^ ) is proportional 
to ff ■ r, and since fi_ lies in the plane of the quarks, 
it is proportional to a linear combination of Yii{r) and 
li_i(f), with the junction position r defined relative to 
a (body) z-axis perpendicular to the quark plane. When 
the quarks are in their lowest energy Lq — state in the 
hybrid confining potential, it can be shown that the to- 
tal angular momentum projection with respect to fixed 
axes is ±1 as for the body axes, so that the total or- 
bital angular momentum must be at least L = I. This 
L = 1 configuration is expected to be dominant in the H 
hybrids. 

The junction Hamiltonian in Eq. ( |0.3| ) is invariant un- 
der the inversion of coordinates, -~^ — r; which implies 
li — li, and r — r, so that the fiux wavefunction must 
be a state of good parity. Since ?7_ is a vector in the plane 
of the quarks it can be written as a linear combination of 
the li, with coefficients which are functions of the parity- 
invariant lengths li. It follows that is parity- odd , and 
that the hybrid-baryon wavefunction from Eq. (^^) has 
even parity. 

Permutations Pij of the quark labels exchange the 
quark positions r; <-> Vj, and so 1^ <-> Ij. Since the junc- 



tion Hamiltonian in Eq. (0.3) is symmetric under such 



a relabeling, then Hfiux and Pij must commute. This 
implies that "if and Pij"^ , where 5" is an eigenfunction 
of iJfiux with energy y(ri,r2,r3), must be degenerate. 
Since the baryon B and hybrid H have different ener- 
gies, Pij'^ must be a multiple of VE". As the only one- 
dimensional representations of the permutation group are 
either totally symmetric (S) or totally antisymmetric (A), 
it follows that for the hybrid ^E" is either S (hybrids H^) 
or A (hybrids H^) under quark label exchange. With the 
above flux Hamiltonian, the logical choice is an exchange- 
symmetric baryon trial wavefunction ^ b- 

The color structure of baryons and hybrid baryons is 
motivated by the strong-coupling limit of lattice QCD, 
where the quarks are triplet sources of color joined by 
strings connected to a totally antisymmetric junction, so 
that the color wavefunctions of both kinds of state are 
totally antisymmetric under quark-label exchange. The 
Pauli principle therefore requires that the ground state 
symmetric hybrids have combined flavor-spin wave- 
functions for the quarks which are totally symmetric, 
since these will have symmetric quark orbital wavefunc- 
tions. This implies that flavor-symmetric A states must 
have 5 = |, and that flavor-mixed-symmetry N states 
have S = ^. The hybrids require totally antisymmet- 
ric quark flavor-spin wavefunctions, which are impossible 
for A states, and possible only for N states with S" = i. 
When these flavor-spin configurations are put together 
with the expected dominant L^ ~ 1+ ® 0+ combined 
fiux-quark orbital wavefunction for the H hybrids, the 
result is the set of configurations shown in Table ||. Al- 
though both the present model and the bag model |^ 
predict the presence of seven low-lying hybrid baryons. 
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importance. 



Hybrid Baryon 


L 


S 


(iV,A) 


2S+1 jP 










1 


1 3 

2 ' 2 


iV21 + 
2 


N^3 + 
Jv 2 , 


^ 2 ' 


^ 2 ' 


A*r 




1 


1 
2 


iV21 + 
^' 2 


Ar23 + 

2 








bag model 


1 


1 3 

2 ' 2 


^ 2 ' 


^ 2 ' 


Ar4i + 


iV43 + 






1 


1 
2 


iV2i + 

2 


Ar23 + 

^' 2 









TABLE II. Quantum numbers of low-lying hybrid baryons 
for the adiabatic surface H, degenerate in the absence of 
spin-dependent forces. 
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only the states N'^^'^ and N'^^'^ have the same flavor, 
quark spin S, total angular momentum, and parity as 
low-lying hybrids in the bag model Flavor assign- 
ments for the hybrids are reversed in the bag model, 
due to the different (mixed) exchange symmetry of the 
quark color wavefunctions required by the addition of a 
single colored gluon. These differences emphasize the col- 
lective nature of the gluonic excitation in the flux-tube 
model. With the addition of spin-dependent forces, the 
lowest-lying hybrid states are predicted to be two pairs 
of quark-spin- i nucleon states in the Pn and P13 partial 
waves, with more massive quark-spin-| A states. The 
central conclusion of this letter is that the lightest of 
these states are at 1870±100 MeV, considerably higher 
than previous mass estimates from the bag model and 
QCD sum rules [| of about 1500 MeV. 

If hybrid baryons obey suppression of certain decay 
modes similar to the decay selection rules for hybrid 
mesons p^ , they may be distinguishable from conven- 
tional baryons in the same mass range on the basis of 
their strong decays; in addition, their electromagnetic 
couphngs are expected to be distinctive There are 

baryon states predicted by quark potential models to 
be present in these partial waves at these energies which 
are missing from the analyses. A signal for the presence 
of hybrid baryons would be the discovery, in analyses of 
data expected from the new experiments, of more states 
than predicted by models which constrain the glue to be 
in its ground state. Furthermore, in order to understand 
experiments which are designed to find these missing 
baryons, states with excited glue must be considered. An 
obvious place to look for signals for hybrid baryons would 
be in electro- and photoproduction of p and to in Hall B 
at TJNAF; there are also planned ttA^ scattering exper- 
iments by the Crystal Ball collaboration at the new 
D-line at BNL which will examine the final states Nr], 
Np and Nuj in the 2 GeV mass region. 

The decays —* ppco and — > pprj have been ob- 
served |2^] with branching ratios of roughly 10~^. Since 
gluonic hadron production is expected to be enhanced 
above conventional hadron production in the glue-rich 
decay of the ^p, it is possible that a partial wave analysis 
of the pu or pr] invariant masses would yield evidence 
for hybrid baryons. Future work in this area at BEPC 
and an upgraded r-charm factory could be of critical 
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